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ABSTRACT 

Based on the latest SNe la data provided by Hicken et al. (2009) with 
using MLCS17 light curve fitter, together with the Baryon Acoustic Oscilla- 
tion(BAO) and strong gravitational lenses(SGL), we investigate the constraints 
on the dark energy equation-of-state parameter w in the flat universe, especially 
for the time-varying case w(z) = Wq + w z z/(l + z). The constraints from SNe 
data alone are found to be: (a) {VL M ,w) = (0.358, —1.09) as the best-fit results; 
(b) (wo, w z ) = (— 0.73^q;97, 0.84^}q 6 3 4 ) for the two parameters in the time-varying 
case after marginalizing the parameter Qm] (c) the likelihood of parameter w z 
has a high non-Gaussian distribution; (d) an extra restriction on f2 m is necessary 
to improve the constraint of the SNe la data on the parameters (wq, w z ). A 
joint analysis of SNe la data and BAO is made to break the degeneracy between 
w and Qm, and leads to the interesting maximum likelihoods w = —0.94 and 
w z = 0. When marginalizing the parameter Qm, the fitting results are found 
to be (wq,w z ) = (— 0.95io!i8j 0-4lia96)- After adding the splitting angle statis- 
tic of SGL data, a consistent constraint is obtained (f2 M ,w) = (0.298,-0.907) 
and the constraints on time-varying dark energy are further improved to be 
(wo,w z ) = (— 0.92^°;^, 0.35^q;54), which indicates that the phantom type models 
are disfavored. 
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1. INTRODUCTION 

Analysis of the distance modulus versus redshift relation of type la supernova (SNe la) 
provides a direct evidence that the universe expansion is accelerating in the last few billion 
years(e.g., Perlmutter et al. 1999, Riess et al. 1998, 2004, 2007; Astier et al. 2006; Wood- 
Vasey et al. 2007; Kowalski et al. 2008; Hicken et al. 2009; Kessler et al. 2009). This 
cosmic image is also supported by many other cosmological observations, like the Cosmic 
Microwave Background(CMB)(Hinshaw et al. 2009; Komatsu et al. 2009), the Baryon 
Acoustic Oscillation (BAO) measurement (Eisenstein et al. 2005; Cole et al. 2005; Huetsi 
2006; Percival et al. 2007) and the weak gravitational lenses (Weinberg and Kamionkowski 
2002; Zhan and Knox 2006). Based on the Friedmann equation, the acceleration can be 
explained through introducing a negative pressure component in the universe, named dark 
energy, which is nearly spatially uniform distribution and contributes about 2/3 critical 
density of universe today. To reveal the property of dark energy, the most of studies, either 
theoretical models or experiment data analysis, are focused on its equation-of-state parameter 
w = p/p. Here we shall utilize the latest SNe la data provided by Hicken et al. (2009) 
with using MLCS17 light curve fitter, together with the splitting angle statistic of strong 
gravitational lenses(SGL; Zhang et al. 2009) and the baryonic acoustic oscillations (BAO; 
Eisenstein et al. 2005) to investigate the constraint for the parameter w of dark energy in 
the flat cosmology, especially for the time varying case. 

By far, all observed data are consistent with the ACDM cosmology, with dark energy 
in the form of a cosmological constant A. However, this model raises theoretical problems 
related to the fine tuned value (see e.g. Padmanabhan 2003). Many other theoretical models, 
like quintessence models and phantom model(Ratra and Peebles 1988; Caldwell, Dave, and 
Steinhardt 1998; Caldwell 2002) , reveal that dark energy might be a dynamical component 
and evolves with time. It is usual to parametrize dark energy as an ideal liquid with its 
equation-of-state(EOS) parameter w(z) = w + w z z/(l + z) (Chevallier et al. 2001; Linder 
2003). Conveniently, it includes the case of a constant EOS with (w = w,w z = 0), and 
the ACDM model (w = —l,w z = 0). Then theoretical models can be classified in a phase 
diagram on the (wo,w z ) plane (see e.g. Barger et al. 2006; Biswas and Wandelt 2009). 
Thus the accurate measurement of the parameters (wo,w z ) is helpful for testing a certain 
theoretical models. The current allowed regions of (wo,w z ) given by the most observations 
or their combinations remain surrounding the crucial point (wq = —l,w z = 0), which is 
the common point in the phase diagram of different classified models. Therefore, the final 
judgment of models can not be made and more careful works are still needed. 

The SNe la has homogeneity and extremely high intrinsic luminosity of peak magnitude 
and thus is widely used to measure the cosmological parameters(e.g. Riess et al. 2004; 
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Wood-Vasey et al. 2007; Kowalski et al. 2008; Hicken et al. 2009; Kessler et al. 2009). With 
given density of dark energy in the universe today p(z = 0), the change of equation-of-state 
parameter w will bring the change of its density p(z) at the redshift z and then the distance 
d(z). Inversely, the measurements of the redshift z of SNe la and corresponding distance 
d(z) can constrain the dark energy. In spite of the high accuracy of SNe la measurement, 
its potential of constraining dark energy is not very strong due to the degeneracy of w and 
Qm- Therefore, a combining analysis with other observations is often useful. 

We shall also use the summary parameters of the baryonic acoustic oscillations as re- 
ported in previous studies(Eisenstein et al. 2005). The large-scale correlation function of a 
large sample of luminous red galaxies has been measured in the Sloan Digital Sky Survey 
and a well-detected acoustic peak was found to provide a standard ruler by which the abso- 
lute distance of z — 0.35 can be determined with 5% accuracy, which is independent of the 
Hubble constant h. This ruler is a VLm prior and can be used to constrain dark energy(e.g. 
Porciani and Madau 2000; Huterer and Ma 2004; Chae, 2007). We are going to show that the 
strong gravitational lensing statistic observation can also provide us a useful probe of dark 
energy of the universe. This is because the dark energy affects, mainly through the comoving 
number density of dark halos described by Press-Schechter theory and the background cos- 
mological line element, the efficiency with which dark-matter concentrations produce strong 
lensing signals. Then by comparing the observed number of lenses with the theoretical ex- 
pected result as a function of image separation and cosmological parameters, it enables us 
to determine the allowed range of the parameter w. The constraint process also depends 
on the density profile of dark halos. Here we will use the two model combined mechanism 
to reproduce the observed curve of lensing probability to the image splitting angle (Sarbu, 
Rusin and Ma 2001; Li and Ostriker 2002; Zhang et al. 2009). The redshift of CMB is above 
1000 and far larger than 1, and there is no other observation to fill up this redshift gap, 
thus we would not adopt the CMB data in the present analysis and limit our study on dark 
energy to the redshift region of z ~ 1, which is the characteristic redshift scale of SNe la, 
BAO and SGL statistic. 

In our recent work(Zhang et al. 2009), we have present the constraint on the dark energy 
from the SGL splitting angle statistic. In this paper, by taking the latest analyzed SNe la 
data(Hicken et al. 2009), the baryonic acoustic oscillations(Eisenstein et al. 2005) and the 
CLASS statistical sample(Browne et al.2003), we shall make a joint analysis to constrain the 
dark energy equation of state parameters w, especially for the time- varying parameterization 
w(z) = w + w z z/{\ + z). We mainly highlight two issues which have not previously been 
illuminated. First, we carefully study based on the latest SNe la data the constraints for the 
dark energy EOS w(z) and the influences of Qm- Second, we investigate the joint analysis 
of SNe la data, BAO and SGL statistic in detail. Our paper is organized as follows: Sect. 2 
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shows the constraint by the latest SNe la data on dark energy. Sect. 3 describes the joint 
analysis of the SNe la, the BAO and the SGL statistic, more stringent constraints on dark 
energy are resulted. The conclusions are presented in the last section. 

2. DARK ENERGY CONSTRAINTS BY THE LATEST SNe la DATA 

As the standard candles of the cosmology, the SNe la is used to study the geometry and 
dynamics of the universe with redshift z < 1.7. In determinations of cosmological parameters 
about the accelerating expand and dark energy, the SNe la remains a key ingredient. In 1998, 
the SNe la measurement provided the first direct evidence for the presence of dark energy 
with the negative pressure. Then many SN la observations have been done and the total 
number of SNe la sample increases quickly. The SN la compilations are often consist of 
high-redshift (z ~ 0.5) data set and low-redshift (z ~ 0.05) sample at the same time (e.g. 
Riess et al. 1998; Perlmutter et al. 1999; Wood-Vasey et al. 2007). When combining several 
independent group's SNe la data sets into one compilation, the consistent analysis method 
of light curves and the selection of supernova are crucial. For a certain sample, the different 
light curve fitter and corresponding different selection of supernova can lead to different 
constraints on the cosmological parameters(e.g. Hicken et al. 2009). 

The fitting results of cosmological parameters from different SN la compilations have 
moderate differences. To obtain the consistent and more powerful constraint, researchers 
have made many efforts to deal with the cross-calibration uncertainities when combining the 
different SNe samples. There is a conventional method to combine several group's SNe la 
compilations, namely, by introducing an extra nuisance parameter in the x 2 statistic of every 
used SNe la sample and marginalizing them over in the fit, all \ 2 statistics of samples can be 
summed into one total statistics (see e.g. Barger et al. 2006). The nuisance parameters are 
considered as analysis-dependent global unknown constants in the distances. Although this 
combined mechanism is wildly adopted, the so-called analysis-dependent unknown constant 
is just an averaged effect of analysis-dependent uncertainties. 

Kowalski et al. (2008) provided the Union data set, a compilation of 307 SNe la dis- 
covered in different surveys. The heterogeneous nature of the data set have been reflected 
and all SNe la sample are analyzed with the same analysis procedure. In the Union data 
set, all SNe la light curve are fitted by using the spectral-template-based fit method of Guy 
et al. (2005) (also known as SALT). There are other light curve fitters used in literatures, 
such as SALT2 (Guy et al. 2007), MLCS2k2 (Jha, Riess, and Kirshner 2007) with R v = 3.1 
(MLCS31) and MLCS2k2 with R v = 1.7 (MLCS17). Hicken et al. (2009) compared these 
light curve fitters and found that SALT produces high-redshift Hubble residuals with sys- 
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tematic trends versus color and larger scatter than MLCS2k2, and MLCS31 overestimates 
host-galaxy extinction while MLCS17 does not. For a certain SNe la, the analysis outcomes 
of different light curve fitters are not equal. Here we choose the SNe la compilation provided 
by using MLCSf 7 light curve fitter with the best cuts Ay < 0.5 and A < 0.7 to constrain 
the dark energy. 

In the flat universe, the Friedmann equation are given by 

h(z)/h = y/n M {i + z y + n DE {z) 

J (f -fi M )(l + ^) 3(1+w) for constant w, 

DE{Z) ~ \ (l-fi M )(l + z)3(i+™o+™*) e -3™W(i+ z ) for w ^ = WQ + Wz _|_ ; W 

with Hubble constant H = fOO h km s _1 Mpc -1 . The influence of cosmological parameter 
w is focused on the dark energy density VLde(z) and then the Luminosity distance d^, which 
is defined as 

Analysis of the distance modulus versus redshift relation of SNe la can give us the 
information about the cosmological parameters. Distance estimates of SNe la are derived 
from the luminosity distance, di = (C/AnJ-') 1 ^ 2 where C and T are the intrinsic luminosity 
and observed flux of the SNe la, respectively. It is usual to introduce the apparent magnitude 
m and absolute magnitude M. From the definition of the distance moduli \x = m — M, we 
have 

y u = 51ogd i /Mpc + 25. (3) 

Using Equations (CQ), (j2J) and (J3j), we can relate the parameter w with the measured redshift 
z and distance moduli fx{z) of SNe la data. Since parameter Hq is irrelevant for the SN only 
data, the likelihood of the SNe la analysis can be determined from a x 2 statistic 



^(n M , ro) = ^ MM^l^Lg ( 4) 



where subscript i denotes the zth SNe la data and <jj is the observed uncertainty. fi° and 
/i T are the observed and theoretical distance moduli, respectively. 

Let us first discuss the constraints for the constant w case. Using the Powell mini- 
mization method(Press, et al 1992), we minimize the likelihood function of the parameters 
(Qm,w) and find that the coordinate of the best-fit point is (Ojw, w) = (0.358,-1.09). 
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Though the results of SNe la data and WMAP observation are consistent in the statistical 
meaning, it is noted that this best-fit value of VIm is large, in comparison with Qm — 0.258 
of the concordance cosmology provided by WMAP five year data(Komatsu et al. 2009). 
Figure [T] shows the likelihoods of parameter Qm and w, in which the maximum likelihood 
points are located at VLm = 0.36 and w = —0.88, respectively. It is interesting to notice that 
the parameter w is restricted to be from —2.0 to —0.5 and Qm is less than 0.5. 

We now focus on the time- varying model w(z) = w + w z z/(l + z). Figure [2] shows 
the contours of (w ,w z ), the best-fit point is (w ,w z ) = (—0.73,0.84) after marginalizing 
the parameter Qm- It is seen that the SNe la data alone have a poor constraint power 
on the parameter w z . In figure [3], we plot the contours of two parameters (Qm,wo) after 
marginalizing w z , the best-fit point is found to be (f2 M ,u7 ) = (0.45,-0.68). It is shown 
that when Qm increases from 0.3 to 0.45, the allowed region of parameter wq is enlarged 
quickly. For a smaller Qm < 0.3, it leads to a much better constraint for the parameter 
wo: wo ~ —1.4 ~ —0.6. Figure 0] gives the contours of (Qm,w z ) after marginalizing wq, the 
best-fit point is (Qm,w z ) = (0.44, —4.63). It is noticed that when VLm increases from 0.34 to 
0.5, the allowed region for the parameter w z is enlarged rapidly. For a smaller VLm < 0.34, 
it also leads to a much better constraint for the parameter w z : w z ~ —3.0 ~ 2.5. From the 
figure [3] and figure HI it indicates that an extra restriction on Qm is necessary to improve the 
constraint of the SNe la data on the parameters wq and w z . Figure [5] shows the likelihoods 
of parameters Wo and w z , in which the maximum likelihood points are located at w = —0.8 
and w z = 0.4, respectively. It can be seen that the parameter w$ is limited in the region 
—2.5 < Wq < 0.5 and the likelihood of parameter w z has a high non-Gaussian distribution. 

3. JOINT ANALYSIS OF SNe la DATA, BAO AND SGL STATISTIC 

As we have shown in the previous chapter, over 300 SNe la observed so far are not 
sufficient for determining the cosmo logical parameters, especially for w and w z . Many 
surveys(e.g. the Dark Energy Survey and Pan-STARRS) are proposed to obtain the SNe la 
sample with enlarged number and improved precision. Here we are going to constrain the 
dark energy through the combination of the SNe la data, the baryon acoustic oscillations as 
well as the SGL splitting angel statistic. 
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3.1. Baryon Acoustic Oscillations 

In the relativistic plasma of the early universe, ionized hydrogens (protons and electrons) 
are coupled with energetic photons by Thomson scattering. The plasma density is uniform 
except for the primordial cosmological perturbations. Driven by high pressure, the plasma 
fluctuations spread outward at over half the speed of light. After about 10 5 years, the 
universe has cooled enough and the protons capture the electrons to form neutral Hydrogen. 
This decouples the photons from the baryons, which dramatically decreases the sound speed 
and effectively ends the sound wave propagation. Because the universe has a significant 
fraction of baryons, these baryon acoustic oscillations leave their imprint on very large scale 
structures (about lOOMpc) of the Universe. 

The measurement of baryon acoustic oscillations was first processed by the Sloan Digital 
Sky Survey (SDSS; York et al. 2000) and Eisenstein et al. (2005) studied the large-scale 
correlation function of its sample, which is composed of 46, 748 luminous red galaxies over 
3816 square degrees and in the redshift range 0.16 to 0.47. The typical redshift of the sample 
is at z = 0.35. The large-scale correlation function is a combination of the correlations 
measured in the radial (redshift space) and the transverse (angular space) direction (Davis 
et al. 2007). Thus, the relevant distance measure is modeled by the so-called dilation scale, 
D v (z) = [D\(z)z/H(z)] 1 ^ 3 , with comoving angular diameter distance D A (z) = f* dz'/H(z'). 
The dimensionless combination A(z) = D v (z)^ VLmH% / z has no dependence on the Hubble 
constant h and is found to be well constrained by the SDSS data at z — 0.35. A standard 
ruler is provided as (Eisenstein et al. 2005) 



M 



[H( Zl )/H }V3 



i r zi 

Zl Jo 



dz 



Zl J H(z)/H 



2/3 



0.469 ±0.017, (5) 



where z\ = 0.35. This ruler is a Qm prior and can be used to constrain dark energy. Then 
the statistic is given by 

2/ \ [A(n M ,w) - 0.469] 2 

x (n M , w) = (6) 



3.2. SGL Splitting Angle Statistic 

The CLASS statistical sample has provided a well-defined statistical sample with N = 
8958 sources. Totally iVj = 13 multiple image gravitational lenses have been discovered and 
all have image separations A6 < 3" (Browne et al. 2003). The SGL statistics are sensitive 
to the equation-of-state parameter w of dark energy, which influences the number density 
of lens galaxies and the distances between the sources and lens. The probability with image 
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separations larger than AO for a source at redshift z s on account of the galaxies distribution 
from the source to the observer can be obtained by (Schneider et al. 1992) 

P(>A0) = / — -(1 + z) 3 n(M,z)a(> AO) dMdz , (7) 
Jo Jo dz 

where M is the mass of a dark halo, Dl is the proper distance from the observer to the lens, 
n(M, z) is the comoving number density of dark halos virialized by redshift z with mass M 
and a is the cross section for two images with a splitting angle > AO. 

According to the Press-Schechter theory, the comoving number density with mass in the 
range (M, M + dM) is given by 

n(M,z)dM = ^f(M,z)dM. (8) 

with the matter density of universe today po — ^MPcrit,o and the critical matter density at 
present p C rit,o = SHq/(8tiG). f(M,z) is the Press-Schechter function, and we shall utilize 
the modified form by Sheth and Tormen (1999) 



0.383 S. dA 



A 2 xO.3- 



0.7075? 



x exp 



0.707 f5c 
A 



(9) 



v^F A 2 dM 

A 2 (M,z) = J —A k (k,z)W 2 {kr) (10) 

where A is the variance of the mass fluctuations(Eisenstein and Hu 1999) and parameter 
5 c (z) is the linear overdensity threshold for a spherical collapse(Wang and Steinhardt 1998; 
Weinberg and Kamionkowski 2002). 

For different density profile of dark halo, the lensing cross section a can be calculated 
out based on the lensing equation. We shall use the combined mechanism of SIS and NFW 
model to explain the whole experimental curve of strong gravitational lensing statistic. For 
that a new model parameter M c was introduced by Li and Ostriker (2002): lenses with mass 
M < M c have the SIS profile, while lenses with mass M > M c have the NFW profile. Then 
the differential probability is given by 

dP/dM = dP S is/dM$(M c — M) + dP NFW /dM${M - M c ) 

where d is the step function, $(x — y) = 1, if x > y and otherwise. As the splitting angle 
A# is directly proportional to the mass M of lens halos, the contribution to large A# of SIS 
profile is depressed by M c . The lens data require a mass threshold M c ~ 1O 13 /i _1 M , which 
is consistent with the halo mass whose cooling time equals the age of the universe today. 
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The likelihood function of the SGL splitting angle statistic is defined as 



L(w) = (l-p(w)) N - N >l[q i 



w). 



(11) 



p(w) and qiiw) represent the model-predicted lensing probabilities and the differential lensing 
probabilities, respectively. They are related to P in Equation (jTJ) by an integration 



p(w) = P obs (> A9) 




and 



q(w) = 



dP ohs (> Afl) 
dA6 




B ip(z s )dzdz s 



p d 2 P{> AO) 

B dAOdz ^s)dzdz s 



(12) 



(13) 



B is the magnification bias and can be found in our previous work(Zhang et al. 2009). (p(z s ) 
is the redshift distribution of sources. Here we take the Gaussian model by directly fitting 
the redshift distribution of the subsample of CLASS statistical sample provided by Marlow 
et al. (2000), which is given by (Zhang et al. 2009) 



exp 



-(z a - a)' 
2A 2 



(14) 



with A, = 1.6125; a = 0.4224: A = 1.3761. 



3.3. Joint Analysis and Numerical Results 

In this section, we will investigate the constraint on the cosmological parameters from 
the joint analysis of (SNe + BAO) and (SNe + BAO + SGL), respectively. For the two(or 
three) independent observations, the likelihood function of a joint analysis is just given by 

L = £sNe X -^BAO ( x -^SGl) 

= exp(-xlNc/2) x exp(-xl AO /2) (xL SG l). (15) 

The statistic significance Xbao an< l XsNe can t> e obtained by using Equations and (JJJ), 
respectively. Lsgl is the likelihood function of SGL statistic and can be obtained by using 
Equation fllip . For SGL data, we shall integrate the parameter h from 0.4 to 0.9. 

Let us first discuss the constraints for the constant w case. In figure we show the 
constraints on Qm and the constant w from the joint analysis of (SNe + BAO + SGL). For a 
comparison, the results of (SNe +BAO) have been shown as dotted lines. The best fit result 
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is (CIm,w) = (0.29, —0.91). The 95% C.L. allowed regions of constant w and Qm are found 
to be: -1.06 < w < -0.77 and 0.25 < Q M < 0.34. Comparing with the results of (SNe + 
BAO) case, it is seen that the results have only slight differences and the fitted w is found 
to be slightly smaller after adding the SGL data. 

Figure [7] plots the likelihoods for the parameters Qm and w from the joint analysis of 
(SNe + BAO) and (SNe + BAO + SGL), respectively. The maximum likelihood points are 
located at Q M = 0.29 and w = -0.88 for (SNe + BAO) and Q M = 0.296 and w = -0.91 for 
(SNe + BAO + SGL). It is interesting to find that the parameters w and Qm are restricted 
to the range: -1.17 < w < -0.67 and 0.23 < Q M < 0.37. 

After marginalizing the parameter Qm, we obtain the constraint on (wq, w z ) in figure 
IE] from the joint analysis of (SNe + BAO) and (SNe + BAO + SGL), respectively. The 
crosshairs mark the best-fit point (w , w z ) = (—0.95, 0.41) for the (SNe + BAO) case and 
(w ,w z ) = (-0.92,0.35) for the (SNe + BAO + SGL) case. For the (SNe + BAO) case, the 
95% C.L. allowed regions for the parameters u>o and w z are found to be: —1.22 < u> < —0.66 
and -0.92 < w z < 1.59. For the (SNe + BAO + SGL) case, the 95% C.L. allowed regions 
for the parameters w and w z are found to be: —1.10 < w < —0.72 and —0.55 < w z < 1.32. 
After adding the SGL data, the constraint on the parameter wo is improved moderately, but 
for the parameter w z , the allowed region decreases by near half. The extra constraint power 
on the time- varying w(z) obtained through adding SGL data is due to the larger redshift 
< z < 3.0 of the galaxies in CLASS observational sample, in comparison with the reshift 
range of SNe data < z < 1.5 and the redshift of BAO z = 0.35. It can be seen for the both 
cases that: (a) the most allowed region of w z is above w z = 0; (b) in comparison with the 
cosmological constant (wo,w z ) = (—1.0,0.0), the joint analysis for both cases favors more 
positive (wq,w z ); (c) in comparison with the results of SNe la data alone, the constraint on 
w z is much improved and Wq also gets better constrained. 

Figure |9] plots the likelihoods of parameters wq and w z from the joint analysis of (SNe 
+ BAO) and (SNe + BAO + SGL), respectively. For (SNe + BAO) case, the maximum 
likelihood points are located at w = —0.94 and w z = 0. Note that w z = implies a constant 
equation-of-state of dark energy. For (SNe + BAO + SGL) case, the maximum likelihood 
points are found to be w = —0.91 and w z = 0.34. We see that the parameters w and w z 
are restricted to be: —1.20 < wq < —0.67 and —1.0 < w z < 2.0. 
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4. CONCLUSIONS 

We have carefully investigated, based on the latest SNe la data, BAO and SGL statistic, 
the constraint on the equation-of-state parameter w of dark energy, especially for the time 
varying cases in the flat universe. The influences of the matter density Qm on the fitting 
results are carefully demonstrated. The typical redshift measured by the three kinds of 
observations is z ~ 1 and far smaller than the redshift of CMB involved, their constraints 
on the parameter w are effective and significant only for the redshift region z < 1.5. 

The influence of the equation-of-state parameter w on the density p(z) of dark en- 
ergy in the universe and the distance d(z) makes SNe la data a powerful probe of dark 
energy. In this paper, we have utilized the latest 324 SNe la data provided by Hicken et 
al. (2009) using MLCS17 light curve fitter with the best cuts A v < 0.5 and A < 0.7 to 
carefully investigate the constraint on the equation-of-state parameter w of dark energy. 
For the constant w case, the best-fit results for the two correlated parameters are found 
to be (Qm,w) = (0.358,-1.09). It is seen that Qm is somewhat large in comparison with 
Qm = 0.26 of the concordance cosmology provided by WMAP five year data(Komatsu 
et al. 2009); note that using a different parameterization of dark energy, an alterna- 
tive analysis (Huang et al. 2009) presented a best-fitted result VL M = 0.446 from SNe la 
data, which is even larger but still consistent with our result at 95% C.L. For the time- 
varying case, after marginalizing the parameter Q M , we have obtained the fitting results 
(w ,w z ) = (— 0.73^o;g7, 0.84lJ 6 g 4 ), which indicates that (a) the SNe la data alone have only 
a poor constraint power on the parameter w z , an extra restriction of Qm is necessary, so 
that the constraint of SNe la on the parameters Wq and w z can be much improved; (b)the 
likelihood of parameter w z has a high non-Gaussian distribution. 

The summary parameter of BAO can provide a standard ruler by which the absolute 
distance of z = 0.35 can be determined with 5% accuracy. This ruler can be a Qm prior and 
has been used to constrain dark energy. The strong gravitational lensing (SGL) statistic is 
a useful probe of dark energy. Through comparing the observed number of lenses with the 
theoretical expected result, it enables us to constrain the parameter w. We have used the 
latest SNe la data together with the BAO (and the CLASS statistical sample) to constraint 
dark energy. For the constant w case, the results obtained from (SNe + BAO) and (SNe + 
BAO + SGL) only have a slight difference. We have shown that: (a) for the (SNe + BAO) 
case, the best fit results of the parameters w) are (0.287, —0.885) and for the (SNe + 
BAO + SGL) case, the best fit point is (Sl M ,w) = (0.298,-0.907); (b) the fitting results 
are found to be Q M = 0.29^;^ and w = -0.91^;^ for the (SNe + BAO + SGL), which 
are consistent with the ACDM at 95% C.L.; (c) the most allowed region of parameter w is 
above the line w — — 1. Comparing with the fitting results from the SNe la data alone, we 
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have found: (a) the allowed region at 95% C.L. for Vt M is reduced to one-fifth; (b) the best 
fit value of w is almost not changed but its variance is reduced very much. 

For the time- varying case w(z) after marginalizing (Qm), we have obtained the fit- 
ting results (w ,w z ) = (-0.951^,0.41^) for the (SNe + BAO) case and (w ,w z ) = 
(-0.9218;^, O.SbtoZ) for the (SNe + BAO + SGL) case. It has been seen that the adding of 
the SGL data makes the constraints on parameter (wo,w z ) to be much improved. For both 
cases, the most allowed region of w z is above w z = 0, which indicates that the data from 
the three observations (SNe + BAO + SGL) disfavor the phantom type models. Comparing 
with the fitting results from the latest SNe la data alone, we have observed that: (a) the best 
fit values for wq are decreased by over 0.2 and the variances are approximately reduced to 
one- fourth; (b) the best fit values of w z are decreased by 0.49 and the variances are reduced 
to one-twelfth. 

In conclusion, the joint analysis of the latest MLCS17 data set given by Hicken et al. 
(2009), summary parameters of BAO and SGL data have provided an interesting constraint 
on the equation-of-state parameter w of dark energy, especially for the time-varying case 
with parameters (w ,w z ). A large number of SNe la samples with reduced systematical 
uncertainties in the near future, together with possible new observations on BAO and SGL 
statistic, would be very useful to understand the properties of dark energy. 
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Fig. 1. — The likelihoods of the parameters Qm and w. The maximum likelihood points are located 
at Q.m = 0.36, and w = —0.88, respectively. 
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Fig. 2. — The (wo,w z ) contours of SNe la data alone. The crosshairs mark the best-fit point 
(w ,w z ) = (-0.73,0.84). 




Fig. 3. — The (Am, wo) contours of SNe la data alone after marginalizing w z . The crosshairs mark 
the best-fit point (Q M ,w ) = (0.45,-0.68). 
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Fig. 4. — The (Am, w z ) contours of SNe la data alone after marginalizing wq. The crosshairs mark 
the best-fit point (£l M ,w z ) = (0.44,-4.63). 
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Fig. 5. — The likelihoods of the parameters wo and w z . The maximum likelihood points are located 
at wq = —0.8 and w z = 0.4, respectively. It can be seen that the likelihood of parameter w z has a 
high non-Gaussian distribution. 
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Fig. 6. — 68% C.L. and 95% C.L. allowed regions of (Qm,w) from the joint analysis of (SNe + 
BAO + SGL) (solid lines) in comparison with the joint analysis of (SNe + BAO) (dotted lines). 
The best-fit result from (SNe+ BAO + SGL) is (Sl M ,w) = (0.29, -0.91). 
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Fig. 7. — The likelihoods of the parameters Q,m and w from the joint analysis of (SNe + BAO) 
and (SNe + BAO + SGL), respectively. The maximum likelihood points are located at Q.m = 0.29 
and w = -0.88for (SNe + BAO) and Q M = 0.296 and w = -0.91 for (SNe + BAO + SGL). 



-23 - 




Fig. 8. — the 68% C.L. and 95% C.L. allowed regions of (wo,w z ) from the joint analysis of 
(SNe + BAO) and (SNe + BAO + SGL), respectively. The crosshairs mark the best-fit point 
(w ,w z ) = (-0.95,0.41) for the (SNe + BAO) case and (w ,w z ) = (-0.92,0.35) for the (SNe + 
BAO + SGL) case. 



-24 - 




-1.4 -1.2 -1.0 -0.8 -0.6 -2-1012 



Fig. 9. — the likelihoods of parameters u>o and w z from the joint analysis of (SNe + BAO) and 
(SNe + BAO + SGL), respectively. For (SNe + BAO) case, the maximum likelihood points are 
located at wo = —0.94 and w z = 0.0, respectively. For (SNe + BAO + SGL) case, the maximum 
likelihood points are wq = —0.91 and w z = 0.34, repectively. 



